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Abstract 

This note is devoted to continuity results of the time derivative of the solution to the one- 
dimensional parabolic obstacle problem with variable coefficients. It applies to the smooth fit 
principle in numerical analysis and in financial mathematics, ft relies on various tools for the 
study of free boundary problems: blow-up method, monotonicity formulae, Liouville's results. 
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1 Introduction 

Consider a parabolic obstacle problem in an open set. We look for local properties, which 
do not depend on the boundary conditions and the initial conditions, but only depend on 
the equation in the interior of the domain. Consider a function u with a one-dimensional 
space variable x in Qi(0) where by Q r (Po) we denote the parabolic box of radius r and 
of centre Pq — (xo,to) : 

Q r (P ) = {(x,t) Gl 2 , \x-x \<r, \t-t \<r 2 } . 

Assume that u is a solution of the one-dimensional parabolic obstacle problem with 
variable coefficients: 

J a(x, t)u xx + b(x, t)u x + c(x, t)u-ut = f(x, t) ■ l {u>0} a.e. in Qi(0) 
1 b>0 a.e. in Q x (0) 
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where u t ,u x ,u xx respectively stand for ^4r, and lr u >o} is the characteristic 

function of the positive set of u. Here the free boundary T is defined by 

r = (3{ii = o})ngi(o) . 

To simplify the presentation, we assume that the coefficients 

a, b, c and / are C 1 in (x, i) , (2) 

but Holder continuous would be sufficient in what follows. 

A natural assumption is that the differential operator is uniformly elliptic, i.e. the 
coefficient a is bounded from below by zero. If we do not make further assumptions on 
a and on /, we cannot expect any good property of the free boundary T. Suppose that: 

3 6 > 0, a(x,t) > 5, f(x,t) > 6 a.e. mQi(O) . (3) 

Up to a reduction of the size of the box (see 0]), any weak solution u of JJJ has a bounded 
first derivative in time and bounded first and second derivatives in space. Assume there- 
fore that this property holds on the initial box: 

|wt(x,t)|, |ita;(a;,£)| and \u xx (x,t)\ are bounded in Qi(0) . (4) 

This problem is a generalisation to the case of an operator with variable coefficients 
of Stefan's problem (case where the parabolic operator is d 2 /dx 2 — d/dt). Stefan's 
problem describes the interface of ice and water (see ^| ED |5]). The problem with 
variable coefficients arises in the pricing of american options in financial mathematics 

(see OHmHUEHElCL]) 

If P is a point such that u(P) > 0, by standard parabolic estimates Ut is continuous 
in a neighbourhood of P. On the other hand if P is in the interior of the region {u = 0}, 
Ut is obviously continuous. The only difficulty is therefore the regularity on the free 
boundary T. By assumption ut is bounded but may be discontinuous on T. The regularity 
of u t is a crucial question to apply the "smooth-fit principle" which amounts to the 
C 1 continuity of the solution at the free boundary. This principle is often assumed, 
especially in the papers dealing with numerical analysis (see P. Dupuis and H. Wang [SJ 
for example). 

In a recent work L. Caffarelli, A. Petrosyan and H. Shagholian 5 prove the C°° 
regularity of the free boundary locally around some points which are energetically char- 
acterised, without any sign assumption neither on u nor on its time derivative. This 
result holds in higher dimension but in the case of constant coefficients. We use tools 
similar to the ones of [S] and the ones the last author developed previously for the elliptic 
obstacle problem in ^2]. Our main result is the following: 
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Theorem 1.1 (Continuity of u t for almost every time) Under assumptions 
for almost every time t, the function u t is continuous on Qi(0). 

This result is new, even in the case of constant coefficients. The continuity of Ut cannot be 
obtained everywhere in t, as shown by the following example. Let u{x,t) = max{0, — t}. 
It satisfies u xx — ut — l{ u >o} an( i its time derivative is obviously discontinuous at t = 0. 
If additionally we assume that ut > we achieve a more precise result: 

Theorem 1.2 (Continuity of u t for all t when ut > 0) Under assumptions 
*/ u t > in Qi(0) then u t is continuous everywhere in Qi(0). 

The assumption that Ut > can be established in some special cases (special initial 
conditions, boundary conditions, and time independent coefficients). See for example 
the results of Friedman [7j, for further results on the one-dimensional parabolic obstacle 
problem with particular initial conditions. 

In Section |3 we introduce blow-up sequences, which are a kind of zoom at a point of 
the free boundary. They converge, up to a sub-sequence, to a solution on the whole space 
of the obstacle problem with constant coefficients. Thanks to a monotonicity formula for 
an energy we prove in Section|21that the blow-up limit is scale-invariant. This allows us to 
classify in Section 0] all possible blow-up limits in a Liouville's theorem. Then we sketch 
the proof of Theorcm ll.il We even classify energetically the points of the free boundary 
into the set of regular and singular points. In Sectional we prove the uniqueness of the 
blow-up limit at singular points. Then we give the sketch of the proof of the Theorem 
11.21 For further details we refer to 0]. 

2 The notion of blow-up 

Given a point Pq = (xo,to) on the free boundary T, we can define the blow-up sequence 

by 

u k (x,t) = uiX0 + £ f 0+£H) >e >0. (5) 

Roughly speaking the action of this rescaling is to zoom on the free boundary at scale e 
(see figure P). 

By assumption, u(Pq) = 0. Because u is non-negative, we also have u x (Pq) = 0. 
Moreover u E Po has a bounded first derivative in time and bounded second derivatives in 
space. For this reason, using Ascoli-Arzela's theorem, we can find a sequence (e n ) n which 
converges to zero such that (u £ p Q ) n converges on every compact set of I 2 = R x x M.t to 
a function u° (called the blow-up limit) and which a priori depends on the choice of the 
sequence (e 
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Figure 1: Blow-up 



The limit function u° satisfies the parabolic obstacle problem with constant coeffi- 
cients on the whole space-time: 

a(P )u° xx - «t = /( P o) • hu>0} m M 2 . 

By the non-degeneracy assumption (3), it is possible to prove that € = 0}. 

To characterise the blow-up limit u , we need to come back to the original equation 
satisfied by u and to obtain additional estimates. In order to simplify the presentation 
we make a much stronger assumption on u: assume that u is a solution on the whole 
space-time of the equation with constant coefficients o = 1, f = 1, b = and c = 0: 

u xx -u t = 1{„> } in K 2 . (6) 

Without this assumption, all tools have to be localised. See 0] for more details. 

3 A monotonicity formula for energy 

For every time t < 0, we define the quantity 

£(*;«) = j ^— t (\u x (x 1 t)\ 2 + 2u(x,t))~^u 2 (x,t)^G(x,t)dx 
where G satisfies the backward heat equation G xx + Gt — in {t < 0} and is given by 
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Theorem 3.1 (Monotonicity formula for energy) Assume thatu is a solution of^Bjj. 
The function £ is non- increasing in time for t < 0, and satisfies 

j t £{t- u) = -— L3 J \Cu(x, i)| 2 G(x, t) dx (7) 

where 

Cu(x, t) = —2 u(x, t) + x • u x {x, t) + 2t ■ u t (x, t) . 
A similar but different energy is introduced in jS] . 

Corollary 3.2 (Homogeneity of the blow-up limit) Any blow-up limit u° of(u E p o ) n 

defined in Jty). satisfies 

u a (\x, X 2 t) = A 2 u°(x, t) for every x e E, t < 0, A > . (8) 
Proof. We prove it in the case Po = 0. The crucial property is the scale-invariance of £: 

£{e 2 n t-u) =£(t;u £ Q ") . 

Taking the limit e n 0, we get 

£(Q~;u):= lim £(t; u) = £(t; u°) for every t<0. (9) 

r^O r<0 

From the monotonicity formula J7J), we get Cu (x,t) = for t < 0. This implies the 
homogeneity JSJl of u° for t < 0. □ 

4 A Liouville's theorem and consequences 

Let 

v+{x,t) = -(max{0,x}) 2 , 
V-(x,t) = — (max{0, — x}) 2 , 

and for m £ [—1,0] 

1 + m n 



v m (x,t) = < 



m t + 



tV m [^) >0 if 0<t<C m -x 2 , 



if i > C m • x 2 , 



where the coefficient C m is an increasing function of m, satisfying C m = if m = — 1, 
and C m = +oo, if to = 0. The precise expression of V m is given in [3]. In particular we 
get z;_i(x,i) = max{0, — t] and Uo(^>i) = \x 2 - 
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Theorem 4.1 (Classification of global homogeneous solutions in E 2 ) Letu ^ 
be a non-negative solution of satisfying the homogeneity condition J^. Then u° is 
one of v+, w_ or v m for some m G [—1,0]. 




u+ V- i) m ,mg(— 1,0) V-i Vq 

Figure 2: Solutions of Theorem 14. II 

Similar versions of this theorem are also proved in [5]. 

Theorem II .21 is a consequence of Theorem 14.11 Every blow-up limit satisfies w° < 0. 
A more detailed analysis leads to 

liminf u t (P) < . (10) 

P^Po 

From the assumption u t > we so infer that u t = 0. 

We also have an energy criterion to characterise points of the free boundary 

Theorem 4.2 (Regular and singular points) Letu be a solution o/@). Then either 
£(0-;u) =V2 or £(Q-;u) = y/2/2. 

In the first case (i.e. £(0~;u) = \/2) Pq is called a singular point. Otherwise (i.e. 
£(0~; u) = V2/2), Pq is a regular point. 

Proof. By © we have £ (0~; u) = £ (— 1; u°). Blow-up limits have been classified in Theo- 
rem l4.ll A simple calculation gives £ (— 1; v + ) = £(—l;v + ) = V2/2, and £(—l; v m ) = \/2 
for every m£ [-1,0]. □ 

5 A monotonicity formula for singular points 

One of the crucial idea of ^2] can be adapted to the parabolic framework. 

Theorem 5.1 (Monotonicity formula for singular points) Let u be a solution of 0) 
and assume that Pq — is a singular point. For any m G [—1,0] the function 

t ^ <P Vm (t;u) = [ i (u(x,t) - v m (x,i)f G(x,t) dx , t<0 (11) 

is non-increasing. 
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As a consequence lim T ^ , r<o ® Vm u) := $ , ' m (0 ; u) is well denned. 

Corollary 5.2 (Blow-up limit at singular points) Under the assumption of Theo- 
rem \5.1\ there exists a real m G [—1,0] such that any blow-up limit of u at is equal to 
v m . 

Proof. Consider two sequences (e n ) n and (e„)„ converging to such that (uQ n ) n and 
( u o")n respectively converge to two blow-up limits u° — v m , for some m £ [0, 1], and vP. 
Using the scale-invariance we get 

$^(-l;<«)=^(-4t;n) and *»™ (-1; = *«»(-^t;u) . 

Passing to the limit in the scale-invariance we obtain 

= *«™(_l;ti°) = * o ™(0-;tt) = *°"*(-l;fi ) . 

This proves that it = u m = u°, i.e. the uniqueness of the blow-up limit. □ 

Heuristically the singular points whose blow-up limit is v rn with m G [—1,0) have a 
free boundary with horizontal tangent in the (cc, f)-plane: this almost never occurs. On 
the other hand the singular points with vq as blow-up limit and the regular points have 
a blow-up limit which satisfies = 0. This last argument can be refined to show that 
the time derivative of u is continuous at such points, and consequently for almost every 
time. This proves Theorem ll.il 
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